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1.1 Introduction 

In the twenty five years that have passed since Berry published his landmark paper |1], the geometric 
phase has turned out to be a very unifying concept in physics. Various phase shifts and rotation 
angles both in classical and quantum physics have been proven to originate from the geometry of the 
underlying parameter space. One of the first examples was given by Pancharatnam [2 ] who discovered 
that the phase shift due to a cyclic transformation of the polarization of an optical field is equal to half 
the enclosed area on the Poincare sphere for polarization states. Other optical examples of geometric 
phases are the phase shift that arises from the variation of the direction of the wave vector of an optical 
field through a fiber and the phase that is associated with the cyclic manipulation of a squeezed 
state of light [4]. The Gouy phase shift, which is due to the variation of the beam parameters (the 
beam width and the radius of curvature of the wave front) of a Gaussian optical beam, can also be 
interpreted geometrically Q. 

In analogy with the geometric phase for polarization (or spin) states of light, van Enk has proposed 
a geometric phase that arises from cyclic mode transformations of paraxial optical beams carrying 
orbital angular momentum [6]. The special case of isotropic first-order modes is equivalent to the 
polarization case Q and, as was experimentally demonstrated by Galvez et. al., the geometric phase 
shift acquired by a first-order mode that is transformed along a closed trajectory on the corresponding 
Poincare sphere also equals half the enclosed surface on this sphere JSJ. Similar experiments were 
performed with second-order modes [9], in particular to show that exchange of orbital angular mo- 
mentum is necessary for a non-trivial geometric phase to occur iflQl . However, in the general case of 
isotropic modes of order N, the connection with the geometry of the N + I dimensional mode space 
is not all obvious. 

In this chapter, we present a complete and general analysis of the phase shift of transverse optical 
modes of arbitrary order when propagating through a paraxial optical set-up, thereby resolving this 
issue. Paraxial optical modes with different transverse mode indices (n, m) are connected by bosonic 
ladder operators in the spirit of the algebraic description of the quantum-mechanical harmonic oscilla- 
tor and complete sets of transverse modes \u nm ) can thus be obtained from two pairs of ladder operators 
ifTTIl . We show that the geometries of the subspaces of modes with fixed transverse mode numbers n 
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and m, which are closed under mode transformations, are all carbon copies of the geometry under- 
lying the ladder operators. We fully characterize this geometry including both the generalized beam 
parameters, which characterize the astigmatism and orientation of the intensity and phase patterns of 
a Gaussian fundamental mode, and the degrees of freedom associated with the nature and orientation 
of the set of higher-order modes. We find a dynamical and a geometric contribution to the phase shift 
of a mode under propagation through an optical set-up, which both have a clear significance in terms 
of this parameter space. 

The material in this chapter is organized as follows. In the next section we briefly summarize 
the operator description of paraxial way optics. We discuss its group-theoretical structure, which is 
essential for our ladder-operator approach, and show how paraxial ray optics emerges from it. In 
section 11.31 we discuss how complete basis sets of transverse modes can be obtained from two pairs 
of bosonic ladder operators. We discuss the transformation properties of the ladder operators, and, 
thereby, of the modes and characterize the ten degrees of freedom that are associated with the choice 
of a basis of transverse modes. Two of those degrees of freedom relate to overall phase factors of the 
ladder operators and, therefore, of the modes. In section [L4l we show that the variation of these phases 
under propagation through a set-up originates from the variation of the other parameters. We discuss 
an analogy with the Aharonov-Bohm effect in quantum mechanics and show that both contributions to 
the phase shift are geometric in that they are fully determined by the trajectory through the parameter 
space. However, only the geometric contribution relates to the geometry of this space. Section [T31 is 
devoted to the specific, but experimentally relevant, case of mode transformations of non-astigmatic 
modes. In the final section, we summarize our results and draw our conclusions. 



1.2 Canonical description of paraxial optics 

1 .2.1 Position and propagation direction as conjugate variables 

A monochromatic paraxial beam of light that propagates along the z direction is conveniently de- 
scribed by the complex scalar profile u{p, z), which characterizes the spatial structure of the field be- 
yond the structure of the carrier wave exp(ikz - icot). The two-dimensional vector p = (x, y) T denotes 
the transverse coordinates. The electric and magnetic fields of the beam can be expressed as 

E(p,z,t) = Rz{E eu(p,z)e ikz - ioJt ) (1-1) 

and 

B(p, z, = Re (e- x e)u(p, zV fe -^} , (1.2) 

where Eq is the amplitude of the field, e is the transverse polarization, e z is the unit vector along the 
propagation direction and co = ck is the optical frequency with c the speed of light. The slowly varying 
amplitude u{p, z) obeys the paraxial wave equation 

hl + 2ik^-)u(p,z) = 0, (1.3) 

where = d 2 /dx 2 + d 2 /dy 2 is the transverse Laplacian. Under the assumption that the transverse 
variation of the field appears on a much larger length scale than the wavelength, this description of 
paraxial wave optics is consistent with Maxwell's equations in free space |[T2l . 
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The paraxial wave equation (11.31) has the form the Schrodinger equation for a free particle in two 
dimensions. The longitudinal coordinate z plays the role of time while the transverse coordinates 
p(x,y) constitute the two-dimensional space. This analogy allows us to adopt the Dirac notation of 
quantum mechanics to describe the evolution of a classical wave field ifTBll . In the Schrodinger picture, 
we introduce state vectors \u(z)) in the Hilbert space L 2 of square-integrable transverse states of the 
wave field, where the z coordinate parameterizes the trajectory along which the field propagates. The 
states are properly normalized (u(z)\u(z)) = 1 for all z and the field profile in real space can be ex- 
pressed as u{p,z) - (p\u(z)). Just as in quantum mechanics, the transverse coordinates may be viewed 
as a hermitian vector operator p = (x, y) J acting on the Hilbert space. The derivatives with respect 
to these coordinates constitute canonically conjugate operators. Rather than the conjugate transverse 
momentum operator -id/dp, which has the significance of the normalized transverse momentum of 
the field, it is convenient to construct the propagation-direction operator by dividing the transverse 
momentum operator by the normalized longitudinal momentum k. Thus, we obtain the hermitian vec- 
tor operator 9 = (& x ,&y) T = -(i/k)(d/dx,d/dy) T . The transverse position and propagation-direction 
operators obey the canonical commutation rules 



where the indices a and b run over the x and y components. In analogy with quantum mechanics, we 
introduce the transverse field profile in propagation-direction representation 



which is the two-dimensional Fourier transform of u(p, z) and characterizes the transverse propagation- 
direction distribution of the field. 

In geometric optics, a ray of light is fully characterized in a transverse plane z by its transverse 
position p and propagation direction 9, which are usually combined in the four-dimensional ray vector 
& T = (p T ,# T ). The operator description of paraxial wave optics may be viewed as a formally quan- 
tized (wavized) description of light rays, where p and 9 have been replaced by hermitian operators 
p and 9 that obey canonical commutation rules (11.41 ) and l/k = X plays the role of fi lfl4l . These 
operators are conveniently combined in the ray operator £• = (p , 9 ). In analogy with quantum me- 
chanics, where the expectation values of the position and momentum operators have a clear classical 
significance in the limit H — > 0, a paraxial wave field reduces to a ray in the limit of geometric optics 
X — > 0. Its transverse position and propagation direction in the transverse plane z are characterized by 
the expectation values (u(z)\p\u(z)) and (u(z)\9\u{z)). 

1.2.2 Group-theoretical structure of paraxial wave and ray optics 

Both the diffraction of a paraxial beam under free propagation, as described by the paraxial wave 
equation (I1.3I) . and the transformations due to lossless optical elements can be expressed as unitary 
transformations |w ut) = U\u m ) on the transverse state of the field. In general, a unitary operator can 
be expressed as 



[pa,k9 b ] = iS ab , 



(1.4) 




(1-5) 



U{{aj\) = e 



(1-6) 



3 



S. J. M. Habraken and G. Nienhuis 



where {a/} is a set of real parameters and {tj} a set of hermitian generators, i.e. fj = Tj. In the 
present case of paraxial propagation and paraxial (first-order) optical elements, the generators are 
quadratic forms in the transverse position and propagation-direction operators. This is exemplified by 
the paraxial wave equation d 1 -3b . which in operator notation takes the following form 

|-Mz)> = -h 2 \u{z)) (1.7) 
oz 2 

and is formally solved by 

/ ik Z 9 2 \ 

Mz)> = expl —J WO)). (1.8) 

This shows that that free propagation of a paraxial field is generated by k9 2 /2, which is obviously 
quadratic in the canonical operators. Since the ray operator £ has four components, the number of 
squares of the operators is four while the number of mixed products is \Z\ = 6, which gives a total of 
ten quadratic forms. They are hermitian and can be chosen as 

Ti = x 2 , T 2 = f , r 3 = xy , T A = k - (M x + & x x) , T 5 = - (y& y + & y y) , 

T 6 = kx&y , T 7 = k0 x , T 8 - k 2 $ x $ y , T 9 = k 2 & 2 x and T lQ = k 2 b 2 . (1.9) 

In terms of these generators, free propagation of a paraxial beam (11.81 ) is described by 

l(f 9 + f 10 )z 



\u(z)) = exp 



2k 



MO)) . (1.10) 



The mixed product fg appears in the generator of free propagation through an anisotropic medium, 
i.e. a medium in which the refractive index depends on the propagation direction 9. In that case 
the propagator can be expressed as exp(-/^0 T N~ 1 §z/2), where N is a real and symmetric matrix that 
characterizes the quadratic variation of the refractive index with the propagation direction. If the 
anisotropy of the refractive index is not aligned along the & x and & y directions, this transformation 
involves fg. A thin lens astigmatic lens imposes a Gaussian phase profile. The unitary transformation 
that describes it can be expressed as 

/ ikp T F~ l p\ 

l«out> = exp k n > , (1.11) 



where F is a real and symmetric 2x2 matrix whose eigenvalues correspond to the focal lengths of the 
lens while the corresponding, mutually perpendicular, eigenvectors fix its orientation in the transverse 
plane. In the general case of an astigmatic lens that is not aligned along the x and y directions, this 
transformation involves the generators t\, f 2 and T3. A rotation of the beam profile in the transverse 
plane can be represented by 

\u mt ) = e^-^u) , (1.12) 

where 7g - f-j - -i(xd/dy - yd/dx) is the orbital angular momentum operator and <f> is the rotation 
angle. The operators T4 and f$ generate transformations that rescale a field profile along the x and y 
directions respectively, i.e. 

u oat (x,y,z) = <p\u oat (z)) - (pk nogfc)f4+nog(c ' )f5 | Min (z)> - u in (c x x,c,y,z) . (1.13) 
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Physically speaking, such transformations correspond to the deformation of a field profile due to 
refraction at the interface between two dielectrics with different refractive indices. 

From the canonical commutation relations (11.41) . it follows that the commutator of any two gen- 
erators (11.91 ) is a linear combination of the generators. In mathematical terms, the algebra of the 
generators is closed, which means that [f^, t{\ = i £ m gki m T m with real structure constants gki m . We 
shall prove that the unitary transformations (11.61 ) with the generators (11.91 ) form a ten-parameter Lie 
group. For reasons that will become clear this group is called the metaplectic group Mp(4). 

Since the states \u(z)) are normalized, the expectation values (u(z)\p\u(z)) and (u{z)\0\u(z)) have the 
significance of the average transverse position and the average propagation direction of the field. A 
special property of the unitary transformations in equation (11.61) with the quadratic generators given 
by (11.91 ), is that the Heisenberg transformation of the vector operator £ T = (p T ,8 T ) is linear, 

so that it can be expressed as 

u\{aj})lQ({aj}) = M({aj})b , (1.14) 

where M({aj}) is the 4x4 ray matrix that describes the transformation of a ray £ T = (p T ,6 T ) under 
the optical element that is described by the state-space operator U({aj}). The denning properties of 
the position and momentum operators, i.e. that they are hermitian and obey canonical commutation 
rules (11.41) . are preserved under this unitary Heisenberg transformation. It follows that M({aj}) is real 
and obeys the identity 

M T ({aj})GM({cij}) = G with G = ( ^ J J, (1.15) 

where and 1 denote the 2x2 zero and unit matrices respectively, so that G is a 4 x 4 matrix. This 
identity (II. 14b ensures that the operator expectation values (u(z)\$-\u(z)) of the transverse position and 
propagation direction transform as a ray, i.e. trace out the path of a ray when the field propagates 
through an optical set-up. This shows how paraxial ray optics emerges from paraxial wave optics and, 
as such, the identity (11.141 ) may be viewed as an optical analogue of the Ehrenfest theorem in quantum 
mechanics iTTSl . The manifold of rays *■ constitutes a phase space in the mathematical sense. The real 
and linear transformations on this manifold that obey the relation (11.151 ). or, equivalently, preserve the 
canonical commutation rules (11.41 ). are ray matrices. The product of two ray matrices is again a ray 
matrix so that ray matrices form a group. The group of real 4x4 ray matrices, which preserve the 
bilinear form * T GJ , where * and J are ray vectors, is called the symplectic group S p(4, R). The 
term symplectic, which is a syllable-by-syllable translation of the Latin 'complex' to Ancient Greek 
and literally means 'braided together', refers to the fact that a phase space is a joint space of position 
and propagation direction (momentum). The 4x4 ray matrices in S p(4, R) emerge from a set of 
unitary state-space transformations, which, as one may prove from equation (11.141) . constitute a group 
under operator multiplication. As was mentioned already, this group is called the metaplectic group 
Mp(4). For real rays &,J€ K 4 , the products z^Gi, and J T GJ vanish. The product & T GJ does not 
vanish and is obviously conserved under paraxial propagation and optical elements. It is called the 
Lagrange invariant lTT6l [T71 and has the significance of the phase-space extent of a pair of rays z, and 
^ . Conservation of this quantity is an optical analogue of Liouville theorem in statistical mechanics. 
The commutators of the quadratic generators Tj and the position and propagation-direction oper- 
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ators are linear in these operators, so that we can write 

- i[fj, l\ = Jjl , (1.16) 

where the 4x4 matrices Jj are real. Explicit expressions of these matrices are given in appendix 
ll.AI Applying equation (11.141 ) to infinitesimal transformations immediately shows that the ray matrix 
corresponding to the unitary state-space operator in equation (11.61 ) is given by 

M({aj}) = e-^jVjJj . (1.17) 

Equation (11.161 ) provides a general relationship between the generators {fj} of the unitary state-space 
transformations (11.61 ) and the generators {Jj} of the corresponding ray matrices (I1.17I ). By applying 

equation (11.151 ) to infinitesimal transformations, one finds that the generators obey J J -G + GJ; = 0. 

j J 

Moreover, from equation (11.161 ) one may prove that 

[[t i ,tj],l] = [J i ,Jj]f. . (1.18) 
Using the Lie algebra [T k , T t ] - i £,„ gu m T m we find that [J k , //] = - £ 

m SkimJm- This proves that the 

metaplectic and symplectic groups are homomoiphic, i.e. for every U e Mp(4) there is a correspond- 
ing M e S p(4, R). The reverse of this statement is not true; a ray matrix M fixes a corresponding 
transformation U up to an overall phase. The homomorphism is an isomorphism up to this phase. 

By using equation (11.151 ) and the expressions of the unitary transformations (11.10I ). (11-1 lb . (11.121 ) 
and (11.131 ) or, equivalently, the relation between (11.161 ) the sets of generators {fj} and {Jj} and the 
definition (11.171 ) of the ray matrices, one finds the 4x4 ray matrices that describe propagation, a 
thin lens, a rotation in the transverse plane and the rescaling of a beam profile due to refraction at 
the interface between two dielectrics. These ray matrices, some of which have been given explicitly 
in references lfT8l [T9l , generalize the well-known ABCD matrices to the case of two independent 
transverse degrees of freedom GUI . 

The group-theoretical structure that we have discussed in this section can easily be generalized to 
the case of D spatial dimensions. In that case there are 2D canonical operators. These give rise to 2D + 
(^2) - 2D 2 + D linearly independent quadratic forms, which generate state-space transformations that 
constitute the metaplectic group Mp(2D). The corresponding ray matrices obey the 2D-dimensional 
generalization of equation (11.151 ) and form the corresponding symplectic group S p(2D, R). In case 
of a single transverse dimension, the three hermitian quadratic forms can be chosen as x 2 , k{x& x + 
& x x)/2 and k 2 d- 2 x . In the analogous case of the quantum-mechanical description of a particle in three 
dimensions, the number of quadratic forms is twenty-one. 

1.3 Basis sets of paraxial modes 
1.3.1 Ladder operators 

As a result of the quadratic nature of the generators (11.9I ). a, possibly astigmatic, Gaussian beam pro- 
file at the z - input plane of a paraxial optical will retain its Gaussian shape in all other transverse 
planes z of the set-up. This is the general structure of a transverse fundamental mode. Complete sets 
of higher-order transverse modes that preserve their general shape under paraxial propagation and 
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paraxial optical elements can be obtained by repeated application of bosonic raising operators a p (0) 
in the z = plane ETTl . In the present case of two transverse dimensions, we need two independent 
raising operators so that p - 1,2. Both the raising operators and the corresponding lowering operators 
a p (0) are linear in the transverse position and propagation-direction operators p and B. Their transfor- 
mation property under unitary transformations e Mp(4) follows from the requirement that acting with 
a transformed ladder operator on a transformed state must be equivalent to transforming the raised or 
lowered state, i.e. 

dlWout) = Ct#l"in> = t>S ( >in> - (I- 19) 

In view of the unitarity of U, this requires that 

€ut = . (1.20) 

Since the generators dl.91 ) are quadratic in the position and propagation-direction operators, the ladder 
operators preserve their general structure and remain linear in these operators under this transforma- 
tion (1 1 -20b - Moreover, their bosonic nature is preserved so that they obey the commutation rules 

[a p (z),al(z)] = 5 pq (1.21) 

in all transverse planes z of the optical set-up if (and only if) they obey bosonic commutation rules 
in the z = plane. When the fundamental Gaussian mode |woo(z)) is chosen such that the lowering 
operators give zero when acting upon it, i.e. ai(z)|woo(z)) = «2(z)l"oo(z)> = 0, the commutation rules 
(11.211) guarantee that the modes 

\u nm (z)) = -pL= (a\(z)) n (al(z)) m |woo(z)> , (1-22) 
ynlml 

form a complete set in all transverse planes z. For a given optical system, the complete set of modes is 
thus fully characterized by the choice of the two bosonic ladder operators a p (0) in the reference plane 
z = 0. 

In reference [18 ], we have shown that, in the special case of an astigmatic two mirror-cavity, the 
ladder operators, and thereby the cavity modes, can be directly obtained as the eigenvectors of the 
ray matrix for one round trip inside the cavity. In the present case of an open system, we are free to 
choose the parameters that specify the ladder operators in the z = input plane. A convenient way to 
do this is to choose an arbitrary ray matrix Mq € S p(4, R). This ray matrix can be chosen independent 
of the properties of the optical system, and of the ray matrices that describe the transformations of 
its elements. However, as we shall see, a necessary and sufficient restriction is that Mq has four 
eigenvectors p for which the matrix element p^Gp does not vanish. It is obvious that this matrix 
element is purely imaginary so that the eigenvectors must be complex. Since Mq is real, this implies 
that for each eigenvector p p also p* p is one of the eigenvectors so that the eigenvectors come in two 
complex conjugate pairs, obeying the eigenvalue relations Mop p = A p p p and Mop* p = A* p p* p , with 
p = 1,2. Without loss of generality we can assume that the matrix elements p P Gp p are positive 
imaginary. Then we can write 

p\Gp p = 2i and p P Gp p = , (1.23) 
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where p = 1,2. The first relation can be assured by proper normalization of the eigenvectors, whereas 
the second follows from the antisymmetry of G. By taking matrix elements of the symplectic identity 
M^GM = G, we find the relations 

A* p A q /j.lG/i q = /UpGfiq and A p A q fi p Gfj. q = n p Gfi q . (1.24) 

Assuming that the two eigenvalues A\ and A2 are different, we conclude that 

n\Gn2 = and H\Gn 2 = . (1.25) 



When the eigenvalues are degenerate, i.e. A\ = A2, one can find infinitely many pairs of linearly 
independent vectors pi and pj that obey these symplectic orthonormality properties. Following the 
approach discussed in reference lPT8l . we now specify the ladder operators in the z = input plane by 
the expressions 

a P (0) = ^i P Gi and fij,(0) = fy p Gl . (1.26) 

The symplectic orthonormality properties (11.231 ) and (11.25b of the eigenvectors p p and p* p ensure that 
the ladder operators in the input plane obey bosonic commutation relations (I1.21I ). From the general 
transformation property of the ladder operators (11.201 ). combined with the Ehrenfest relation (11.141) 
between U and M, one may show that the ladder operators in other transverse planes z are given by 
same expressions (11.261 ) when p p is replaced by p p (z) = M(z)p p . Here, M(z) is the ray matrix that 
describes the transformation of ray from the z = input plane to the transverse plane z. It can be 
constructed by multiplying the ray matrices that describe the optical elements of which the set-up 
consists and free propagation between them in proper order. The fact that the properties (11.231) and 
(11.251) are conserved under symplectic transformations e S p(4, R) confirms that the ladder operators 
remain bosonic in all transverse planes of the set-up. 

Since the modes are fully characterized by the choice of two complex vectors p p , we expect that 
the expectation values of physically relevant operators can be expressed in terms of these vectors. The 
average transverse position and momentum of the beam trace out the path of a ray. This implies that 
the expectation values (u nm \p\u nm ) and (u nm \B\u nm ) vanish. In appendix I l.B I we prove, however, that 
the expectation values of the generators Tj are, in general, different from zero and can be expressed 
as 

(u nm \fj\u nm ) - j ||« + ^jp\GJjPi + {ni + ^jyU^G/^j ■ (1-27) 

This result generalizes the expression of the orbital angular momentum in twisted cavity modes that 
we derived in reference f22l . It is noteworthy that these properties of the modes are fully character- 
ized by the generators Jj and the complex ray vectors p p , which both have a clear geometric-optical 
significance. 

Finally, it is worthwhile to notice that the results of this section remain valid when the number 
of (transverse) dimensions is different. In particular, the same method gives explicit expressions for 
complete orthogonal sets of time-dependent wave functions that solve the Schrodinger equation of a 
free particle in three-dimensional space. 
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1 .3.2 Degrees of freedom in fixing a set of modes 

We have shown that there is a one-to-one correspondence between the defining properties of ray 
matrix, i.e. that it is real and obeys the identity (1 1 - 15b - and the properties (11.231 ) and (11.251 ) of the 
complex eigenvectors p p that ensure that the ladder operators (11.261 ) are bosonic. This implies that 
all different basis sets of complex vectors p p that obey these identities must be related by symplectic 
transformations, i.e. each of these sets can be written as {Mp p \ U {Mp* p \, with M e S p(4,R) and 
{pp) U {p* p \ the set of complex eigenvectors of a specific ray matrix Mq e Sp(4, R). Since {Mp p ) U 
{Mp* p ) constitutes the set of eigenvectors of MM$M~ X , it follows that the freedom in choosing a set 
of complex vectors that generate two pairs of bosonic ladder operators (11.261 ) is equivalent to the 
freedom of choosing a ray matrix M e S p(4, R). As a result, the number of independent parameters 
associated with this choice is equal to the number of generators of S p(4, R), which is ten. In order to 
give a physical interpretation of these degrees of freedom, we follow the characterization discussed 
in reference 11231 and decompose the complex ray vectors into two-dimensional subvectors so that 
y"p(z) T = ( r p(^)' m terms of these subvectors, the ladder operators take the following form 

a p {z) = >/f (rj(z)0 - t T p (z)p) and a\{z) = ^{^6 - t^p) , (1.28) 

where p = 1,2. An explicit expression of the Gaussian fundamental mode can be given if we combine 
the two-dimensional column vectors r p and t p into 

R(z) = (n(z),r 2 (z)) and T(z) = (fi(z),f 2 (z)) • (1-29) 

The objects R and T take the form of 2 x 2 matrices, but since r p and t p are transverse vectors, R and 
T do not transform as such under ray-space transformations € S p(4, R) nor under transformations on 
the transverse plane. The symplectic orfhonormality properties (11.231) and (11.251 ) of the vectors p p can 
be expressed as 

R f (z)T(z) - T t (z)R(z) = 2i1 and R T (z)T(z) - T T (z)R(z) = , (1.30) 
and hold for all values of z. Now, the fundamental transverse mode in plane z can be written as 

" 00(P ' Z) = V^elR(i) eXP l — ) ' (U1) 

where S = -z'TR -1 . As opposed to R and T, S is a 2 x 2 matrix in the transverse plane and transforms 
accordingly. It can be checked directly that acting upon |moo(z)) with the lowering operators a\(z) and 
«2(z) gives zero. The fundamental mode (11.311) is properly normalized and has been constructed such 
that it solves the paraxial wave equation (11.31) under free propagation. Moreover, one may check that 
it transforms properly under the transformations of optical elements. The second relation in equation 
(11.301) guarantees that S is symmetric. This is obvious when we multiply the relation from the left 
with (R T )~ , and from the right with R -1 . The real and imaginary parts S r and S; of S respectively 
characterize the astigmatism of the intensity and phase patterns. The real part can be written as 
S r = ( - z'TR -1 + i(R t )~ 1 T t )/2. With the first relation (fODl this shows that RS r R 1 - 1 . This leads to 
the identity 

RR 1 - S^ 1 , (1.32) 
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which shows that S r is positive definite. As a result, the curves of constant intensity in the transverse 
plane are ellipses. Moreover, the fundamental mode is square-integrable. Depending on the sign 
of det Si(z) the curves of constant phase in the transverse plane are ellipses, hyperbolas or parallel 
straight lines. Under free propagation, S is a slowly varying smooth function of z. Optical elements, 
on the other hand, may instantaneously modify the astigmatism. The astigmatism of both the intensity 
and the phase patterns is characterized by two widths in mutually perpendicular directions and one 
angle that specifies the orientation of the curves of constant intensity or phase. The total number of 
degrees of freedom that specify the astigmatism, and, thereby, the matrix S, is thus equal to six. 

Two of the remaining four degrees of freedom are related to the nature and orientation of the 
higher-order mode patterns. From equation (11.321 ). we find that R can be expressed as S r _1/ V T , where 

— 1/2 

o~ is a unitary 2x2 matrix. Notice that S r is real and positive so that S r is well-defined. It is 
illuminating to rewrite the complex ray vectors m and /i2 as 

(«*>-(;)-u?ro w .Mo T a )<**). **> 

where ju x = (1,0, i, 0) T and p. y = (0, 1,0, z') T are the complex ray vectors that correspond to the ladder 
operators that generate the stationary states of an isotropic harmonic oscillator in two dimensions. 
The first matrix in the second right-hand-side of this expression (11.331) is the ray matrix that describes 
the transformation of a thin astigmatic lens. It imposes the elliptical or hyperbolic wave front of the 
optical modes on the harmonic oscillator functions. The second matrix has the form of the ray matrix 
that describes the deformation of mode due to refraction. It rescales the modes along two mutually 
perpendicular transverse directions and accounts for the astigmatism of the intensity patterns. The 
third ray matrix involves the complex matrix o~ and obeys the generalization of equation (11.151 ) to 
complex matrices. Since it is complex, however, it is not a ray matrix e S p(4, R). In order to clarify 
its significance, we rewrite equation (11.331) in terms of the ladder operators, which are conveniently 
combined in the vector operator (fii,fl2) T - By using the definition if the ladder operators (11.261 ) and 
the Ehrenfest relation (11.14b . the transformation in equation (11.331) can be expressed as 

I ) - a/| RT8 - tT « " -. A /f-e X p(-^)(sy 2 p + i S;'«8)exp(^) . (1.34) 

The linear combination of the position and momentum operators between the brackets takes the form 
of the lowering operators for an isotropic harmonic oscillator in two dimensions. Again, the 2x2 
matrix S r accounts for the astigmatism of the intensity patterns by rescaling the ladder operators 
and, therefore, the modes they generate. The exponential terms take the form of the mode-space 
transformation for a thin astigmatic lens and impose the curved wave fronts. From right to left, 
the lowering operators (11.341) as well as the corresponding raising operators, first remove the curved 
wave front, then modify the mode patterns and eventually restore the wave front again. The 2x2 
matrix cr is a unitary transformation in the space of the lowering operators a\ and &2 and transforms 
accordingly. It arises from the U (2) symmetry of the isotropic harmonic oscillator in two dimensions 
and accounts for the fact that any, properly normalized, linear combination of the bosonic lowering 
operators yields another bosonic lowering operator. Up to overall phases, to which we come in a 
moment, this transformation can be parameterized as a\ —> r]\a\ + 772^2 and 0,2 — * + il*2® 2 w * tn 
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|?7i| 2 + |?72| 2 = 1. The two obvious degrees of freedom that are associated with the spinor 77 = (771, r]2) T 
are the relative amplitude and the relative phase of its components. Analogous to the Poincare sphere 
for polarization states (or the Bloch sphere for spin- 1/2 states), they can be mapped onto a sphere. 
For reasons that will become clear, this sphere is called the Hermite-Laguerre sphere ifTTT . Since 771 
and 772 are spinor components in a linear rather than a circular basis, this mapping takes the following 
form 



n = 





1 ' 




= V2\ 



e 2 cos j + e 2 sin| 

■ !£ A _;£ ■ A 

-ie 2 cos j + te 2 sin | 



(1.35) 



where & and tp are the polar and azimuthal angles on the sphere. The mapping is such that the north 
pole (# = 0) corresponds to ladder operators that generate astigmatic Laguerre-Gaussian modes with 
positive helicity. The south pole (■& = n) corresponds to Laguerre-Gaussian modes with the opposite 
helicity while the equator (1? = n/2) corresponds to Hermite-Gaussian modes. Other values of the 
polar angle & correspond to generalized Gaussian modes Il24l . The azimuth angle ip determines the 
transverse orientation of the higher-order mode patterns. Since paraxial optical modes are invariant 
under rotations over n in the transverse plane, the mapping in equation (11.351 ) is such that a rotation 
over ip on the sphere corresponds to a rotation of the mode pattern over <p = ip/2. 

The unitary matrix that describes the ladder operator transformation corresponding to the spinor 
77 is constructed as 

m m 

* * 

-12 v l 

where the second row is fixed up to a phase factor by the requirement that cr must be unitary. With this 
convention, the two rows of sigma correspond to anti-podal points on the Hermite-Laguerre sphere. 
Completely fixing the matrix cr e U(2), however, requires 4 independent degrees of freedom. The 
remaining two, which are not incorporated in 77, are phase factors. Any matrix cr e U(2) can be written 
as 

' «ta 



(1.36) 



& 1 

The phase factors exp(/y p ) correspond to overall phase factors in the vectors p. p and, therefore, in the 
ladder operators (1 1 -26b . The vectors fi p can be written as 

H P = e^v p (S,Tf) , (1.38) 

where p — 1, 2 and v p (S, 77) is completely determined by S and 77 according equation (11.331) . cr being 
replaced by o" (?7). Although the vectors v\ and V2 obey symplectic orthonormality conditions (11.231) 
and are, therefore, not independent, the phases \\ and ^2 are independent. From equation (11.371) and 

— 1/2 T 

the fact that R = S r cr it is clear that the argument of det(R) is equal to^i + xi so that the overall 
phase of the fundamental mode (11.311) is given by -ix\ +X2)/2. The overall phases of the two raising 
operators are respectively -\\ and -\2, so that the phase factors in the higher order modes \u nm {z)) 
are given by exp(-7^„ m ) with 

Xnm = \n + -Ui +\m + -)x2 ■ (1-39) 



In a single transverse plane, such overall phase factors do not modify the physical properties of the 
mode pattern. The evolution of these phase under propagation and optical elements, however, can be 
measured interferometrically. 
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The astigmatism of the modes, as characterized by the 2x2 matrix S, can clearly be modified in 
any desired way by appropriate combinations of the optical elements that we have discussed in section 
11.21 As will be discussed in section 11.51 the degrees of freedom associated with the spinor n can be 
manipulated by mode convertors and image rotators. Although we shall see that variation of the phase 
factors exp(ix p ) is, in general, unavoidable when the other parameters are modified, we show here 
that it is possible to construct a ray matrix € S p(4, R) that solely changes these phase factors. Such a 
ray matrix is defined by the requirement that 

M x {\x P )) (/u lf i 2l i\ l i* 2 ) = ( efri/ij 6^1*2 e~ ix 'fi\ e~^fi* 2 ) , (1.40) 

so that the vectors p p and p* p are eigenvectors of M x . The corresponding eigenvalues are unitary. In 
terms of R and T this relation can be expressed as 





R* \ 




R* \ 




) 


( T 


T* )' 


( t 


J* ) 


(o 


C* J 



M x ({Xp))\ t „ 

where 



e 1 ^ 
e* 2 

By using that 

R R* r 1 1 ( -V Rt 



(1.42) 



T T* ) 2i{ T T -R T 
which follows directly from the identities in equation (11.301 ). we find that M x can be expressed as 



(1.43) 



M x {\x P )) = - 



\_( R R*||C \( -V Rt \ _ ]_( -RCTt + R*C*T T RCR^ - R*C*R T \ 
2i\ T T* j( C* j[ T T -R T )~ 2i\ -TCT 1 " + T*C*T T TCR t -T*C*R T J 

(1.44) 

This ray matrix adds overall phase factors exp(±;^- p ) to the eigenvectors fi p and fi* p . It is real and one 
may check that it obeys the identity (11.151 ) so that it is a physical ray matrix € S p(4, R). 

In this section, we have argued that the number of degrees of freedom associated with the choice 
of two pairs of ladder operators that generate a basis set of modes in a transverse plane z is equal to 
number of generators of S p(4, R), which is ten. We have shown that six of those are related to the 
astigmatism of the modes as characterized by a the complex and symmetric 2x2 matrix S. Two of the 
other four are angles on the Hermite-Laguerre sphere that characterize a spinor rj, which determines 
the nature and orientation of the higher-order modes. The remaining two are overall phase of the 
ladder operators. All these degrees of freedom can be manipulated in any desired way by choosing a 
suitable ray matrix e S p(4, R). 



1.3.3 Gouy phase 

In the limiting case of non-astigmatic modes that propagate through an isotropic optical system the 
2x2 matrix S is a symmetric matrix with two degenerate eigenvalues so that it can be considered a 
scalar s = s r +isi. If we choose ctq = 1, the higher-order modes are Hermite-Gaussian. In that case, the 
complex ray vectors are given by pi = (r, 0, t, 0) T and p 2 = (0, r, 0, t) T , with r,t e C. The symplectic 
normalization condition (11.231) implies that r*t — t*r = 2i. The real part s T of s = -it/r determines 
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the beam width w = y/2/(ks r ) of the fundamental mode while the imaginary part Sj fixes the radius 
of curvature of its wave fronts according to R = l/s{. Under free propagation over a distance z, the 
vectors fii and ji2 transform according to 



r + zt 

t 




and // 2 (z) 





r + zt 




(1-45) 



The parameters r, t and s remain scalar and free propagation does not introduce an overall phase 
difference between fi\ and H2 so that 77, or, equivalently o"o, is independent of z. Without loss of 
generality we can choose z = to coincide with the focal plane of the mode, which implies that sel 
so that r*t = —t*r = i. Since s t , and, therefore, R = o-QS r cannot pick up a phase, we find that 



X(z) ~X(P) = arg(— = arctan|yj = arctan|^-j , 



(1-46) 



where zr = ir/t is the Rayleigh range. This is the well-known Gouy phase for a Gaussian mode GUI . 
Since the vectors /u\ and /12 pick an overall phase x(z), the raising operators pick up a phase ~x(z)- 
The phase shift of the higher-order modes (11.221) is then given by exp(-/(n + m + \)x) and depends on 
the total mode number N = n+m only. As a result of this degeneracy, the same expression holds in the 
non-astigmatic case with cr + 1 . In that case, it is still true that the components of 77 are independent 
of z. 

Generalization to astigmatic modes is straightforward only if the modes have simple astigmatism 
and if the orientation of the higher-order mode patterns is aligned along the astigmatism of the fun- 
damental mode. In that case, the vectors /u p pick up different Gouy phases and the components of 77 
are independent of z. As will be discussed in section 11.51 this is not true in the case of non-astigmatic 
modes that propagate through an optical set-up with simple astigmatism. In the more general case 
of modes with general astigmatism that propagate through an arbitrary set-up of paraxial optical el- 
ements, the z dependence of S depends on 77 and vice versa ifTTTl . In this case no simple expressions 
of the Gouy phases can be derived. The phase in equation (11.39b may be viewed as the ultimate 
generalization of the Gouy phase within paraxial wave optics. 

1.4 The geometric interpretation of variation of the phases x nm 
1 .4.1 Evolution of the phases Xnm 



In this section we show that variation of the phase differences Xp between /u p and v p (I1.38I) is, in 
general, unavoidable under (a sequence of) mode transformations that modify the degrees of freedom 
associated with S and 77. From the discussion in the previous section it is clear that the generalized 
Gouy phases were defined such that they vary only under transformations that involve free propaga- 
tion. However, for later purposes, it is convenient to formulate the description of mode transformations 
that give rise to phase shifts in a slightly more general way. 

Suppose that the unitary state-space transformation that describes (a part of) a trajectory through 
the parameter space is given by U(0 - exp(-if^), where f is a (linear combination of the) genera- 
tors) defined in equation d 1 -91 ) and ^ is a real parameter that parameterizes the trajectory. In this case, 
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the f dependent ladder operators (11.201 ) obey the anti-Heisenberg equation of motion 

. (1.47) 

In terms of the complex ray vectors /i p (£) and the ray matrix M(£) = exp(-/£) that corresponds to 
U{£) according to relation (1 1 - 14b . this equation of motion takes the form of a symplectic Schrodinger 
equation and can be expressed as 

dC 

Substitution of yu p (£) = exp(ix p )v p (() yields after dividing by exp(/y p ) 



= -Jfi p (0 ■ (1-48) 



■ dx P m dv P 



i^rVpiO + ^7 = -Jv p (0 ■ (1-49) 



By multiplying from the left with v p G, using the normalization condition v p Gv p = 2i and rearranging 
the terms we find that 

— ■ - - { v„GJv„ + v' G-— 

d £ 2 \ p p p di; 



= -WGJv p + viG-^ ■ (1-50) 



The generator / represents a conserved quantity. Hence, the first term between the curly brackets does 
not depend on the parameter f and the above equation (11.501 ) can be integrated to obtain 

XpiO = \ {{vlGJvp) £ + JT dfvjG^j . (1.51) 

The first term between the curly brackets constitutes a dynamical contribution to the phase shift and 
arises from the fact that / corresponds to a constant of motion. The second term, on the other hand, 
relates to the geometry of the complex ray space and is the natural generalization of Berry's geometric 
phase to this case. In the next section, we derive an equivalent expression from which the geometric 
significance of the phase shifts (11.511 ) is more obvious. 



1.4.2 Analogy with the Aharonov-Bohm effect 

In quantum mechanics, it is well-known that the coupling of a particle with charge q to the magnetic 
vector potential A(r) gives rise to a measurable phase shift (q/ti) J c A ■ dr of the wave function when 
the particle moves along a trajectory C = r(t). This effect occurs even when the magnetic field 
B = V x A vanishes everywhere along the trajectory and is known as the Aharonov-Bohm effect |[25l . 

The physical properties that are associated with the wave function that describes a particle in quan- 
tum mechanics are not affected by the transformation ^(r, t) —> exp(/0(r))tA(r, t). The Schrodinger 
equation is obviously not invariant under this local U{\) gauge transformation. When gauge invari- 
ance of the Schrodinger equation is imposed, the vector potential A(r) arises as the corresponding 
gauge field. In this picture, the Aharonov-Bohm phase is due to the coupling to a gauge field, the 
conserved charge q being the coupling constant. As such it is a direct consequence of the U{\) gauge 
invariance of quantum electrodynamics. 

In this section, we point out an analogy between the generalized Gouy phase and the Aharonov- 
Bohm effect. This gives some new insights in the nature and origin of this geometric phase and 
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allows us to derive an expression of the phase (11.511) from which its origin in the geometry of the 
underlying parameter space is obvious. It is convenient to combine the parameters that characterize 
the eight degrees of freedom that are associated with the matrix S and the spinor rj into a vector 
K = (7?i,7?2> • ■•) • The corresponding differential operator, which is a vector in the eight-dimensional 
parameter space, is defined as = (0/0H\,0/0H2, ■ ■■) ■ 

The physical properties, for example those in equation (11.27I ). of the transverse mode fields (11.22I ). 
which are generated by the ladder operators constructed from the vectors p p , are not affected by 
transformations of the type 

H P -» M% p , (1.52) 

where p = 1,2. This symmetry property can be thought of as local U(l) ® U(l) gauge invariance. 
The ray matrix e S p(4, R) that describes such gauge transformations (11.521) figures in equation (1 1.44 1 ). 
As shown in appendix ll.CI the two corresponding real generators J x can be constructed from the 
eigenvectors p p . The vector p\ is an eigenvector of J Xl with eigenvalue -i. Since J Xl is real, the 
complex conjugate vector p* is an eigenvector of J Xp with eigenvalue i. Moreover, J Xl p2 = Jx\^*i ~ 0- 
Similarly, p2 and p* 2 are eigenvectors of J X2 with eigenvalues -i and i and J Xl p\ - J X2 p\ = 0. Since 
invariance under the gauge transformation (11.521) is a local and continuous symmetry, it gives rise 
to conserved Noether charges. The gauge transformations are generated by two different generators, 
hence there are two Noether charges, which can be expressed as v P GJ Xp v p /2 = 1, where the factor 
1/2 arises from the fact that a symplectic vector space is a joint space of position and momentum and 
we have used that J x v p = -i and v\Gv p = 2i. In appendix ll.CI we prove that the corresponding 
state-space generators t Xp can be expressed as (a P a p + a p a^/2 so that the charges of a mode (11.221 ) 
are given by (u nm \f Xl \u nm ) - (n + l/2)and {u nm \t X2 \u nm ) = (m+1/2). Since the gauge transformation 
in equation (11.441) is constructed from the eigenvectors p p , it varies throughout the parameters space. 
As a result, the generators t x can be constructed only locally and vary through the parameter space 
according to the ladder-operator transformation in equation (1 1 -20b . However, since the modes also 
vary, it follows that Noether charges (n + 1/2) and (m + 1/2) of the modes \u nm ) are globally conserved. 
In terms of H and V^, the equations of motion of the vectors p p (11.481 ) can be rewritten as 

oi\ 

These equations are obviously not invariant under the gauge transformations (11.521) . Imposing gauge 
invariance yields the modified equations of motion 

oi\ 

((^ + iA p )v p ). — = -Jv p , (1.54) 

where the gauge fields A p are vector fields in the parameter space of % that are defined by their 
transformation property under the gauge transformations (11.521 ) 

A p ^A p - V^ p . (1.55) 

With these transformation properties, the equation of motion (11.541 ) is manifestly invariant under the 
gauge transformations (1 1 -52b - The general solution of this equation (11.54b can be expressed as 

v,^/ 1 '^, (1-56) 
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where C is a trajectory and [i p solves the equation of motion without the gauge field (I1.53I ). In 
full analogy with the Aharonov-Bohm effect, this shows that the phase difference between fi p and v p 
is due to the fact that the latter is coupled to the gauge field X p . Since we have defined the vectors n p 
so as to include the appropriate geometric-phase factor while they are not coupled to the gauge fields, 
the coupling of v p to the gauge fields removes the geometric phase rather than introducing it. The 
geometric origin of the phases is evident in that they are determined only by the trajectory C and not 
on the velocity dH/d£. By using equation (11.541 ) they can be expressed as 




(1.57) 



which is in obvious agreement with equation (11.51I ). 

In analogy with the Aharonov-Bohm effect, the Noether charges v p GJ x v p /2 = 1 determine the 
strength of the coupling of the vectors fi p to the gauge fields A p . This is consistent with the fact that 
the vectors v p pick up phases Xp- The Noether charges of the modes d 1.221 ). however, are equal to 
n + 1/2 and m + 1/2 and depend on the mode numbers n and m. As a result, the modes \u nm ) couple 
differently to the (corresponding state-space) gauge fields and, therefore, experience different phase 
shifts. This is in obvious agreement with equation dl.391 ). 

The Noether currents (v p GJ x Vp/2)oK/o£ = &R/d£ are uniform throughout the parameters space 

-* 

of K. It follows that the 'physical' fields F a p = d a (A p )p - dp(A p ) a , where the indices a and ft run over 

-* -> 

the parameter-space vector components, are constant so that the gauge fields A p CR) cannot possess 
any non-trivial dynamics. Attributing the generalized Gouy phases Xp to coupling to gauge fields 
X p , which do not have any dynamical properties in their own rights, may seem a bit tautological. 
On the other hand, the analysis discussed here shows that the structure that underlies the generalized 
Gouy phase shifts dl.391 ) is that of a gauge theory. In this picture, the appearance of phase shifts 
under propagation through an optical set-up is the unavoidable consequence of the U(l)xU(l) gauge 
invariance of paraxial wave optics, or, equivalently, of the fact that the mode charges n + 1/2 and 
m + 1/2 are conserved under state-space transformations e Mp(4, R). 

The connection between the gauge invariance as discussed here does not depend on the specific 
structure of the symplectic vector space. Our results as well as the Aharonov-Bohm effect suggest that 
there is a more general connection between local gauge invariance and the appearance of geometric 
phases, see, for instance, reference EoTl . 

1.5 Geometric phases for non-astigmatic modes 
1 .5.1 Ray matrices on the Hermite-Laguerre sphere 

A particulary interesting limiting case of the geometric phases that we discuss in this chapter, are the 
phase shifts due to mode conversions on the Hermite-Laguerre sphere, each point on which character- 
izes a basis set of higher-order modes. We focus on non-astigmatic modes in their focal plane so that 
S can be considered a real scalar s eR. We shall construct ray matrices and corresponding state-space 
operators that solely modify the degrees of freedom associated with the Hermite-laguerre sphere and 
study the geometric phases that arise from such transformations. 
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The azimuth angle <p on the Hermite-Laguerre sphere specifies the orientation in the transverse 
plane of the set of higher-order modes. It can be modified by the rotation operator exp(-iip$3), where 
S3 = k(x& y - }$x)/2 = {f(, - fj)/2 is the corresponding generator. The factor 1/2 accounts for the 
fact that a rotation angle ip on the Hermite-Laguerre sphere corresponds to a (f> = tp/2 in the transverse 
plane. For reasons that will become clear, the ray matrix that describes a rotation in a plane parallel to 
an equatorial plane on the Hermite-Laguerre sphere is denoted M3. It takes the following form 



M 3 (<p) 



P(f) 
P(f) 



cos(jj 


- sin 


sin(f) 


cos 















(?) 
(I) 



cos 



cos I 






(f) 

sin (f) > 



sin 



(1.58) 



where £3 = (76 - Ji)/2 is the corresponding ray-space generator and P e SO(2) is a 2 x 2 rotation 
matrix. When this ray matrix acts on an arbitrary pair of complex ray vectors jj p that obey the identities 
(11.231 ) and (11.251 ). the matrix S transforms according to S — » PSP T . In the present case of scalar S, 
this transformation only modifies the orientation of the mode patterns and does not affect S. 

Another class of transformations that solely act upon the Hermite-Laguerre sphere are those that 
describe mode converters. Mode converters consists of a pair of astigmatic or cylindrical lenses ETl . 
The distance between the lenses and their radii of curvature are chosen such that the Gouy phase shift 
introduces a phase difference # between the eigenvectors jui and fi2 of the converter. If the input and 
output plane of the mode converter are chosen such that they respectively coincide with focal planes 
of the incident and outgoing modes and if the modes are matched to the mirrors so that S is scalar 
and equal to 1 in appropriate units determined by the mirrors, the eigenvectors of the mode converter 
are given by p.\ = (1,0, i, 0) T and fa = (0, 1,0, /) T and their complex conjugates. The ray matrix 
that describes the transformation that introduces a phase difference ■& between fa and fa can then be 
constructed as 

cos(f) sin(f) 

cos(f) 
-sin(f) cos(f) 
sin(f) 

where Si = k(J\ - J<i)/4 + (Jg - J 10) /(4k) is the corresponding generator. Again, one may prove 
easily that this transformation does not affect the astigmatic degrees of freedom if S is scalar and 
equal to 1. The corresponding state-space generator is given by si - k(f\ - f%)/4 + (T9 - fio)/(4k) = 
k(x 2 -f + d-l- #2)/4. 

So far, we have constructed two of the three ray matrices that only modify the nature and orienta- 
tion of the higher-order modes. The third corresponds to a mode converter in a basis that is rotated over 
n/4 in the transverse plane, or, equivalently over n/2 in the equatorial plane of the Hermite-Laguerre 
sphere. The ray matrix that describes such a transformation can be obtained as 



Mi(0) 





-sin(f) 


cos(f) 



(L59) 



M 2 (&) = M 3 (7r/4)M 1 (&)M- 1 (n/4) = e 



COS (|) 




-sin(f) 





cos(f) 
-sin(f) 






sin(f) 
cos(f) 




sin(f) 



cos(f) 



, (L60) 
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where S2 = kJ$/2 + J$/(2k) is the corresponding generator. The corresponding state-space generator 
is given by S2 = ktj/2 + fg/(2k) = k(xy + & x d-y)/2. Since M3 and Mi do not affect the astigmatic 
degrees of freedom if S is scalar and equal to 1, it follows that the same is true for M3. 

By using the canonical commutation relations (11.41 ) and the definitions of the generators Si, $2 and 
S3 in terms of the canonical operators, one may easily show that the generators obey an S U (2) algebra 

[$i,h] = ih (1-61) 

and cyclic permutations. The ray-space generators obey 

[Xi,X 2 ]=I 3 (1-62) 

so that the matrices iE\, 1L2 an d '£3 also constitute an SU(2) algebra. Thus we have obtained both 
a metaplectic and symplectic realization of an S U (2) algebra. This proves the well-known fact that 
SU(2) is a subgroup of Mp(4) and, therefore, of S p(4, R). 



1.5.2 Spinor transformations 



Since the generators Si, X2 and X3 constitute an SU(2) algebra, an arbitrary pair of ray vectors fi p on 
the Hermite-Laguerre sphere can be expressed as a linear combination of the eigenvectors of one of 
these generators. In full analogy with section 1 1 .31 where we introduced the components of 77 as the 
coefficients of the expansion of an arbitrary bosonic lowering operator in terms of the two lowering 
operators for a harmonic oscillator in two dimensions, we can write an arbitrary pair of complex ray 
vectors on the Hermite-Laguerre sphere as 



y"i = mfil + #2 



and 



P-2 = -TJifil + 



(1.63) 



where p\\ - (1,0, i, 0) T and p.2 = (1,0, i, 0) T are eigenvectors of Ei with eigenvalues —i and / respec- 
tively, (they are also eigenvectors of /ho = k{J\ + J 2) 1 4 + (J9 + /io)/(4£) with degenerate eigenvalues 
i). Notice that, analogous to the construction in section [T31 and the mapping in equation (11.351 ). the 
components of 77 are spinor components in a linear rather than in a circular basis. The symplectic or- 
thogonality properties (11.231 ) and (11.251 ) require that the expansions in equation (11.631 ) do not involve 
the complex conjugate vectors /j* . Moreover, they ensure normalization of n such that \n\\ 2 + I772I 2 = 1. 

Since the ray matrices M\, M2 and M3, as defined in equations dl.591 ), (11.601 ) and (11.581 ) only 
modify the degrees of freedom associated with 77, these transformations can be expressed in the two- 
dimensional spinor space. In particular, the transformation described by M3 (11.581 ) can be expressed 
as 

m\i cos (f) - sin (f) 

m ) { sin(f) cos(f) 

where 

' 




- e -i<fT3/2 



m 
m 



(1.64) 



T3 



-1 




(1.65) 



is the corresponding generator. Similarly, the transformations (11.591 ) and (11.601 ) of mode converters 
can be rewritten in terms of the spinor components as 



m 
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(1.66) 
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and 



m 



cos(f) isin(f) \f m ) 
/sin(f) cos(f) ){ m 



where the corresponding spinor generators are given by 



n = | J) and i» = (_°, ■> (1-68) 

As a result of the fact that we have the defined the spinor components with respect to the eigenvectors 
of Hi rather than of E3, the spinor generators n, t% and T3 take the form of Pauli matrices in a rotated 
basis. They also form an SU(2) algebra, i.e. [ti,T2] = n-3 and cyclic permutations. This algebra is 
closed and the matrix transformations of rj on the Hermite-Laguerre sphere that are generated by t\ , 
T2 and T3 are analogous to Jones matrices in polarization optics. 

Since the complex vectors jj\ and are eigenvectors of Si with eigenvalues —i and i respectively, 
the spinor corresponding to a point ((p, 9) on the Hermite-Laguerre sphere can be expressed as 

77(0, 9) = (-^V^+^e-^V^ 2 ( M , (1.69) 

where the factor (— i) ' 2 = exp(-/7r/4) is introduced in order to make this identity consistent with 
equation (11.351) . 

1 .5.3 Mode-space transformations 

In a dimensionless notation, the lowering operators corresponding to the complex ray vectors ju p can 
be expressed as 

b x = + A) and b y = yj^(y + W y ) . (1.70) 

A. J. ^ 4- 

The corresponding raising operators b x and bl generate the set of harmonic-oscillator states in two 
dimensions \v nm ) according to equation (11.221 ). the raising operators being replaced by the harmonic- 
oscillator raising operators. This set corresponds to rj = (1, 0) T , which is on the equator of the Hermite- 
Laguerre sphere. The antipodal point 77 = (0, 1) T gives rise to the same set of modes \v mn ), the mode 
indices being interchanged. The modes corresponding to an arbitrary point on the Hermite-Laguerre 
sphere can be expanded as 

1^(77)) = — = Ubl + rf 2 b]) n (- m bt + m b\T |voo> - (1.71) 
ynlml 

By using that [b x , by] = 0, this can be rewritten as 

(1.72) 

which expresses the transformed state \u nm (rj)) as an expansion in two-dimensional harmonic-oscillator 
states of the same order N = n + ni. Conversely, this result shows that the subspaces of modes of fixed 
order ,/V = n+m are closed under transformations (mode conversions) on the Hermite-Laguerre sphere. 
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Figure 1.1: Intensity and false-color phase patterns of the modes that lie in the tp - plane of the 
Poincare sphere for the non-astigmatic first-order modes |«oi)- The north and a south poles (& = 0,n) 
respectively correspond to Laguerre-Gaussian modes with I — n — m = —\ and I — m — n— 1. On both 
poles p - min(«, m) - 0. The modes on the equator (& = n/2) are Hermite-Gaussian while modes for 
intermediate values of & are generalized Gaussian modes. The color coding in the phase patterns is 
such that the color changes in a continuous fashion from red via yellow, green, blue and purple back 
to red again when the phase changes from to In. 
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In general, the subspace of modes of fixed order N is aN + I dimensional subspace of the Hilbert 
space of transverse states of the field. The unitary transformations on this subspace form the group 
SU(N + 1). Only in the special case of first order modes, the most general unitary transformation is 
equivalent to the SU(2) transformation that figures in equation (11.721) . It follows that the subspace of 
first-order modes is isomoiphic to the Hermite-Laguerre sphere for the ladder operators. This sphere, 
which is an analogue of the Poincare sphere for polarization states (7] , as well as the intensity and 
phase structure of some of the modes that lie on it, is plotted in figure 11.11 In the general case of 
./V > 1, S U (2) is a subgroup of the group S U(N + 1) of unitary transformations on the subspace of 
modes of fixed order N. This accounts for the fact that only specific linear combinations of modes of 
order N are modes as well, i.e. are fully characterized by a set of parameters whose transformation 
through an optical system is governed by the ray matrix M(z). In case of N > 1, the transformation 
in equation (11.721 ) gives rise to a Poincare sphere for each combination (n, m) of the transverse mode 
numbers. Since \u nm {r])) and \u mn (r])) correspond to antipodal points on the same sphere, it follows 
that, depending on the parity of N, only (N + 2)/2 (for even AO or (N + l)/2 (for odd AO of these 
spheres are independent, i.e. not related by rotations over n. They are all isomorphic to the Hermite- 
Laguerre sphere for the ladder operators and may be viewed as analogues of the Poincare sphere for 
specific sets of higher-order transverse modes. The two spheres for second-order modes, as well as 
the intensity and phase patterns of some of the modes that lie on them, are plotted in figure [L2l 

The mode-transformation in equation d 1.72b . together with the matrix representation of the spinor 
transformation that we have discussed above, provides a matrix description of beam transformations of 
non-astigmatic optical modes of arbitrary order. It generalizes the description discussed in references 
E8l|29l , which applies to first order modes. 

By inverting the relations in equation (11.701 ) and their hermitian conjugates, the position and 
propagation-direction operators can be expressed in terms of the ladder operators. Using this re- 

4- A\ 4- A, A, 4- A, A, A, 4- 

suit, the state-space generators can be written as s\ = (b x b x - b y b y )/2, $2 = (b x b y + b x bl)/2 and 

A. 4. A. A. A. -J. 

h - {bxby - b x by)l{2i), which is a Schwinger representation of the SU{2) algebra. Here, the SU(2) 
algebra (11.611 ) is ensured by the boson commutation relations (11.211) . This representation provides a 
complete and closed description of the modes and transformations on the Hermite-Laguerre sphere in 
terms of the ladder operators. 

1 .5.4 Geometric phases and the Aharonov-Bohm analogy 

The spinor 77, as defined by equation (I1.35I) . is completely determined by the azimuthal and polar 
angles on the Hermite-Laguerre sphere. The reverse of this statement is not true; choosing a point on 
the Hermite-Laguerre sphere fixes a properly normalized spinor (ptj = 1 up to an overall phase factor 
so that £ = exp(ix)r](<p, with^ e R. In the limiting case of transformations on the Hermite-Laguerre 
sphere, it follows from equation (11.631) . or from the equivalent expansion in terms of lowering opera- 
tors in equation (11.701) . that the two raising operators a] and at pick up equal but opposite phases ~x 
and x respectively. The modes \u nm {ri)) (11.711 ) pick up a phases exp(-/^-„ m ) with 

Xnm = (n- m)x ■ (1.73) 

Such phases do not modify the physical properties of the modes but their variation under (a sequence 
of) transformations on the Hermite-Laguerre sphere can be measured interferometrically. Analogous 
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Figure 1.2: Intensity and false-color phase patterns of the modes that lie in the ip = plane of 
two the Poincare spheres for the non-astigmatic second-order modes \uu) (this page) and |«02) (next 
page). In the figure on this page, the north and south poles (■& = 0,n) respectively correspond to 
Laguerre-Gaussian modes with l-n-m-m-n-Q and p - min(«, m) = while the modes 
on the equator (& = n/2) are Hermite-Gaussian. The intermediate modes are generalized Gaussian 
modes. In the figure on the next page the north and south poles (& = 0, n) respectively correspond to 
a Laguerre-Gaussian mode with I — n — m - —2 and I — m — n — 2. In both cases p = min(n, m) — 1. 
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(continued) Again, the modes on the equator (& = n/2) are Hermite-Gaussian while generalized 
Gaussian modes appear for intermediate values of &. The color coding in the phase patterns is such 
that the color changes in a continuous fashion from red via yellow, green, blue and purple back to 
red again when the phase changes from to 2n. Both spheres are carbon-copies of Hermite-Laguerre 
spheres on which every point characterizes two pairs of bosonic ladder operators. 
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to the discussion in section 11.41 we shall show that the variation x nas a geometric interpretation 
in terms of the Hermite-Laguerre sphere. We consider (a sequence of) state-space transformations 
that only modify the degrees of freedom associated with the nature and orientation of the higher- 
order modes. The evolution of the ladder operators under such transformations is described by the 
anti-Heisenberg equation of motion (11.471 ) when f is replaced by a generator 5/2, which is a linear 
combination of Si, S2 and S3, and £ parameterizes a trajectory on the Hermite-Laguerre sphere. The 
factor 1 /2 in the generator is introduced for notational convenience. In terms of a spinor the equation 
of motion takes the following form 

* = , ( i.74) 

d£ 2 

where r is the spinor generator that corresponds to S. It is a linear combination of t\, T2 and T3. The 
spinor £ picks up the appropriate phase factor. Substitution of £ = exp(ix)n gives 

• dx dn nrj 

lT} ai; ai = ~~T ' (L75) 

Using that rfrj - 1, this result can be rewritten as 

Tr l1] di- — - (L76) 

The generator r represents a constant of motion so that this result can be integrated to yield 

xi.0 = + 



Jo d(> 



This can also be obtained directly from substitution of the complex ray vectors fi p , as defined by 
equation (11.63I ). in the general expression of the geometric phase shift (11.51I ). The first term in equa- 
tion (11.511 ) arises from the fact that r/2 represents a conserved quantity. The second term constitutes 
the well-known geometric phase shift that is experienced by a spinor when it is transported along 
a trajectory on the Hermite-Laguerre sphere. Analogous to the discussion in section 11.41 both con- 
tributions are geometric in that they are fully determined by the trajectory on the Hermite-Laguerre 
sphere but only the second relates to the geometry of the Hermite-Laguerre sphere. It is natural to 
use spherical coordinates e A = {r sin(#) cos(^o), r sin(#) sin(^o), r cos(#)) T to parameterize points on the 
Hermite-Laguerre sphere. For a closed trajectory that consists of geodesies, the first contribution in 
equation (11.771 ) vanishes OOll . Then, the phase shift (11.771 ) can be rewritten as 

x=t £ * v s = 1 r ^ ' di = 1 § c j ^ ■ d<k ' (L78) 

where is the gradient in spherical coordinates and C - 'R(z) is a closed trajectory on the Hermite- 
Laguerre sphere. In the Aharonov-Bohm picture, this phase shift is due to the coupling of n to a gauge 
field that arises from the U{\) gauge invariance of the spinor dynamics. Comparison with equation 
(11.561) shows that this gauge field is given by 

/=*7 + (f^). (1-79) 
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By using the gradient in spherical coordinates and equation (11.351 ) A can be written as 

-> i cot ■& 

A = — ■ 

The corresponding 'magnetic' field is given by 



(1.80) 



B = ^xA = -^- 2 (1.81) 

and is homogeneous on the Hermite-Laguerre sphere. It may be viewed as the field of monopole 
that is located at the center of the Hermite-Laguerre sphere. By the virtue of Stokes' theorem, the 
geometric phase can be expressed as 

X = i ^ d& = i j)A-d& = i £ M - dS = \j*dCl = |n, (1.82) 

where S is the enclosed surface on the Hermite-Laguerre sphere and Q. is the solid angle. This re- 
sults establishes the well-known connection between the geometric phase acquired by a spinor that is 
transported along a closed trajectory on the Hermite-Laguerre sphere and the enclosed solid angle Q. 
on the sphere. Since we have defined the phase picked up by the higher-order modes as exp(-/^ nm ) 
with^„ m = (n - m)x, the result in equation (11.821 ) has the opposite sign of the analogous result for the 
standard case in which a spinor with positive helicity picks up a phase shift exp(/y). 

The phase shift the modes \v nm ), as given by equation (11.731 ). depends only on the difference 
between the mode numbers n and m. In the Aharonov-Bohm picture, n — m has the significance of the 
topological charge of a non-astigmatic mode \v nm ) and determines the strength of the coupling to the 
(corresponding state-space) gauge field. For modes with equal mode numbers n - m, the topological 
charge vanishes so that they do not couple to the gauge field and, therefore, do not experience a phase 
shift. The orbital angular momentum in non-astigmatic modes \v nm ) can expressed as in - ra)cos# 
Ifm and is proportional to their topological charge. It follows that in the case of a non-astigmatic 
mode, the exchange of orbital angular momentum between the mode and the set-up through which it 
propagates is necessary for a non-trivial geometric phase to occur |[T0l l9l. 

In this section, we have studied the geometric phase that arises from (cyclic) transformations on 
the Hermite-Laguerre for higher-order modes. We have constructed ray matrices that solely modify 
the nature and orientation of the higher-order modes and derived the corresponding spinor and mode- 
space transformations. In terms of the spinor r\ the phase shift due to a (cyclic) transformation takes 
the familiar form of a geometric phase for a spinor. In experimental realizations, mode converters 
consist of pairs of astigmatic lenses in which the degrees of freedom associated with S are employed 
to achieve mode conversion [27]. As a result, there will be an additional contribution to the phase 
shift of the modes. This can be compensated for by measuring the interference between fields that 
have passed the same sequence of mode converters and image rotators but with different relative 
orientations (8). 

1.6 Concluding remarks 

We have explored the parameter space that is associated with the choice of a complete and orthonormal 
set of paraxial optical modes in the transverse plane. Modes are defined as solutions of the paraxial 
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wave equation (11.31 ) that are fully characterized by a set of mode parameters whose variation through 
a paraxial optical set-up is described by the 4x4 ray matrix M(z), which describes the transformation 
of ray r = (p, 6>) T from the z = input plane of the set-up to the transverse plane z. Complete 
sets of transverse modes can be obtained from two pairs of bosonic ladder operators. The ladder 
operators are fully specified by two complex ray vectors p p with p = 1,2, which characterize the 
mode parameters. Their variation through an optical set-up, and, thereby, the variation of the ladder 
operators, can conveniently be expressed in terms of M(z). We have argued that there is a one-to-one 
correspondence between the algebraic properties of the ladder operators and the defining properties of 
a physical ray matrix e 5 p(4, R), i.e. that it is real and obeys the identity (11-1 5b - It follows that all sets 
of modes can be expressed in terms of two pairs of ladder operators and, moreover, that the freedom 
in choosing a set of modes is equivalent to the choice of an arbitrary ray matrix Mq e 5 p(4, R). Since 
S p(4, R) is a ten-parameter Lie group, the number of free parameters associated with this choice is 
equal to ten. A possible physical characterization of these degrees of freedom involves a symmetric 
2x2 matrix S, which characterizes the astigmatism of the phase and and intensity patterns of the 
fundamental mode, and a spinor n, which specifies the nature and orientation of the higher-order 
modes. The matrix S is fully specified by six parameters while characterization of a 77 requires two 
independent parameters, which can be mapped on a Poincare sphere. The remaining two degrees of 
freedom are overall phases of the ladder operators. They do not modify the physical properties of 
the modes in a given transverse plane z. Their variation through an optical set-up, however, gives 
rise to a generalized Gouy phase shift of the modes, which can be measured interferometrically. We 
have shown that both contributions to the variation of the overall phases through an optical set-up, as 
described by equation (11-5 lb . are geometric in that they are fully determined by the trajectory "R{z) 
and do not depend on the velocity &Rjdz- However, only the second contribution in equation (11.511 ) 
relates to the geometry of the parameter space. In the specific case of a closed trajectory on the 
Hermite-Laguerre for non-astigmatic optical modes, the phase shifts of the two raising operators are 
equal but opposite. In full analogy with the Pancharatnam phase, they are equal to half the enclosed 
surface on the sphere. 

It is noteworthy that the overall phases Xp of the vectors p p are in general only unambiguously 
defined in case of a closed trajectory. In particular, in the propagation-direction representation, the 
astigmatism of the fundamental mode uoq(6, z) is fully specified by the symmetric matrix V = S _1 . 
Analogous to the discussion in section [T31 the remaining degrees of freedom can be characterized by 

— 1/2 T 

a unitary 2x2 matrix v, which is defined such that T = V r v . It follows that v and cr are related 

-1/21/2 / — 1 /2 1 /2\ 

by cr - -iv\l r VS r . In general det I V r VS r ) ^ 1 so that defining cr - xo~o and v - x ^0 such 
that o"o and vo have unit determinants, requires different phase matrices x ± x' ■ The phase shift along 
a closed trajectory, however, does not depend on the phase convention used. In the limiting case of 
transformations of non-astigmatic modes in their focal planes, i.e. when S and V can be considered 
real scalars, the phases are also unambiguously defined along an open trajectory. All results presented 
in this chapter are, of course, independent of the phase convention that is chosen. 

We have shown that the symplectic group of ladder-operator transformations S p(4, R) corresponds 
to the metaplectic group Mp{4) of unitary transformations on the Hilbert space of state vectors \u). 
The metaplectic group constitutes a subgroup of the set of all possible unitary transformations. This 
accounts for the fact that only specific linear combinations of paraxial optical modes are modes as 
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well, i.e. are fully characterized by a set of parameters whose variation through a paraxial optical 
set-up fully described by the ray matrix M(z). Each combination (n, m) of the transverse mode indices 
gives rise to a subspace of the Hilbert space of transverse states of the field, which is closed under 
metaplectic transformations. The geometries of these subspaces are all carbon copies of the geometry 
of the symplectic manifold underlying the ladder operators. In the limiting case of mode transfor- 
mations of non-astigmatic modes, the metaplectic group reduces to S U (2) and all those subspaces 
become Poincare spheres, which are all carbon copies of the Hermite-Laguerre sphere for the ladder 
operators. 

We have pointed out an analogy between the Aharonov-Bohm effect in quantum electrodynamics 
and the generalized Gouy effect in classical wave mechanics. This reveals deep insights in the geo- 
metric origin of the latter. The physical properties of the modes (11.221 ) that are generated by two pairs 
of ladder operators are not affected by the gauge transformation described by equation 

(1 1 -52b . or, equivalently (1 1 -44b - Imposing gauge invariance of the equations of motion (11.471 ) or (1 1 -48b . 
gives rise to two gauge fields A p in the parameter space. Analogous to the Aharonov-Bohm effect, 
the geometric phase shift of the ladder operators through an optical set-up is due to the coupling these 
gauge fields. The raising and corresponding lowering operators have pairwise equal but opposite 
topological charges and experience opposite phase shifts. The topological charges of the modes \u nn ), 
i.e. the Noether charges that arise from the gauge invariance of the description of their propagation 
through an optical set-up, are given by n + 1/2 and m + 1/2 and depend on the mode numbers. As 
a result, the modes \u nm ) couple differently to the gauge fields and experience different phase shifts 
given by equation (11.391) . Notice that the above-mentioned subspaces of modes with transverse mode 
indices n and m are all uniquely characterized by their coupling to the two (state-space) gauge fields. 
In the specific case of transformations on the Hermite-Laguerre sphere for higher-order modes, the 
phase shifts of the two lowering operators are equal but opposite. In that case, the phase shift of the 
modes is given by equation (11.731) . In the Aharonov-Bohm picture, the variation of this phase is due 
to the coupling of the spinor rj to a single gauge field A that arises from the U(l) gauge invariance of 
the spinor dynamics. The topological charge of the modes \v nm ) on the Hermite-Laguerre sphere is 
equal n - m and the 'magnetic' field due to the gauge field is uniform on the Hermite-Laguerre sphere. 
It may be viewed as the field of a monopole located at the center of the sphere. 

Although we have focused on the optical case, the mathematical structure that underlies the ladder- 
operator method and the phase shifts that arise from the geometry underlying the ladder operators are 
more general. The ray space (p, 6) is a phase space in the mathematical sense and the operator de- 
scription of paraxial wave optics that we have discussed in section (11.21 ) may be viewed as a formally 
quantized (wavized) description of rays. Although the interpretation is different, all this is in full 
analogy with the quantization of classical mechanics to obtain quantum mechanics. As a result, the 
methods and results of this chapter can be applied to the quantum-mechanical description of wave 
packets. The only restriction for the ladder-operator approach to apply is that the state-space genera- 
tors (or Hamiltonian in the quantum language) are quadratic in the canonical operators p and 9. The 
methods and result in this chapter have been formulated such that it is evident how they can be gener- 
alized to account for more independent spatial dimensions. In the general case of D dimensions, the 
number of generators of Mp(2D) and S p(2D, R) is equal to 2D 2 + D, D 2 + D of which are associated 
with a D X D symmetric matrix that generalizes S. The remaining D 2 parameters specify a unitary 
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matrix e U{D), which generalizes o~, and corresponds to the choice of D orthonormal D-component 
spinors and D overall phase factors. The variation of the phases under propagation (evolution) have a 
geometric interpretation in terms of the other degrees of freedom. 
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1 .A The ray-space generators Jj 

In this appendix we give explicit expressions of the ray-space generators Jj. They are defined by 
equation (II . 16b and correspond to the state-space generators tj as denned in equation ( 11.9b . They are 
given by 
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(1.83) 



1.B Expectation values of the generators f , 



This appendix is devoted to a proof of equation d 1.27b . which expresses the expectation values {u nm \tj\u nm ) 
of the generators tj in equation (11.9b in terms of the corresponding ray-space generators Jj as defined 
by equation (11.16b . We prove this by mathematical induction. The special cases (mooI^/I m oo) involve 
Gaussian standard integrals and can be proven explicitly. A formal proof by mathematical induction 
thus requires showing that the identity (11.27b holds for modes \u nJr \ m ) and \u nm+ \) if it holds for \u nm ). 
In order to prove this, we notice that 



(M n +\ m \T j\u n +\ m } 



n+\ 



(u nm \a\t ja\\u nm ) 



Using that 



(1.84) 



[ f iM = ^{^Gfjl-^Glfj) = ^-f p G[tj, l\ = i^ii T p GJjt , (1.85) 



this can be rewritten as 
/ 1 



,n + 1 



\(u n 



Tja\-iyj-n x GJjl 



&l\ u nm) ~ (u nm \T j\u nm ) — 



n + 1 




'■^H[GJ j{u nm \ia\\u nm ) . 

(1.86) 
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The analogous result may be derived for |« nm+ i) and proving equation (11.271) thus boils down to 
proving that 

- (^j) ^^ P GJ Mnm\^al\u nm ) = ^pGJjUp - X - {/J p GJjfi p y = ^n T pGJjV P , (1.87) 



where we used that G = -G and that J G - -GJ. This expression can be rewritten as 

{u nm \ha p \u nm ) = i(n + vJ^Vp > (1-88) 

which we also prove by mathematical induction. Again, the special case of |woo) can be checked 
explicitly. In order to prove that it is true for \u n+ i m ) and \u nm+ i), we use that 

=4^i G A-[^i G ^= Vf (L89) 



and find 

(u n +\ m \$>a\\u n+ \ m ) - f - - - j (u nm \a\ i-a\a\\u nm ) - 

|jj-j-yj<w„ m |ai f&[& + /-^^yUijaj|«„ m > = (u nm \la[\u nm ) + i '\J^f 1 i ■ (1-90) 

The analogous result may be derived for \u nm+ \). This completes the proof of equation (11.881) and, 
thereby, of equation (1 1 -27b . 

1.C Mode-space operators corresponding to the Noether charges 

In this appendix we construct both the ray-space and the corresponding state-space generators of the 
U(l) ® U(l) gauge transformations. The ray matrix that describes such transformations is given by 
equation (1 1 -44b . To first order in the phases x\ an d Xi tne matrix C (11.421 ) is given by 

_ ( \+ix\ \ ( 1 1 ( i 0) (0 0) 

C= n , . ■ = L , +AT1 L n + *2 n . ■ (1.91) 



(1.92) 



1 + 1X2 ) { 1 ) { ) { / 
Substitution in equation (11.441) then gives 



M^})-1+ 2 ^ _ v _ x hr t f J + 2 [ 



t * T t * T \ 



where r\t\ = r\ ® t\ etcetera are direct vector products. From M x {{x p }) = exp {-Xp^xp) ~ 1 _ XpJx P > 
we find that 



1 ( r„tt + rttl -r,A - r*rj 



XXf 'If. ,"i I (>■»> 



Z. ^ l p Lp T ipi p i p i p Ipl p 

where p = 1, 2. These generators are 4 x 4 matrices in the ray space. By carefully inspecting the form 
of the direct products and the structure of the generators Jj as given in appendix (ll.AI) we find that 

T x» = ~4 VP 0p 1 p + r P 9p h ~ r P 09 r P ~ r P 99 r P + hP 9 r p + l pP 9 r p ~ hPP h ~ f pPP W ■ 



(1.94) 
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This can be rewritten as 



£)UH-*fl(& &)(-*)}• (L95) 

which equals 

(1.96) 

where, we used that £ Gju^ is scalar so that I- G[i p — (£ G/i p ) = -yi p Gt. 
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